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It has long been known that dissipation is a crucial ingredient in the superradiant amplification 
of wavepackets off rotating objects. We show that, once appropriate dissipation mechanisms are 
included, stars are also prone to superradiance and superradiant instabilities. In particular, ultra¬ 
light dark matter with small interaction cross section with the star material or self-annihilation 
can trigger a superradiant instability. On long timescales, the instability strips the star of most 
of its angular momentum. Whether or not new stationary configurations surrounded by scalar 
condensates exist, remains to be seen. 
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I. INTRODUCTION 

Rotating bodies with internal degrees of freedom - 
where energy can be dumped into - display superradi¬ 
ance. This means that the scattering of low-frequency 
waves extracts rotational energy away from the body and 
is used to amplify incoming radiation [H-Q- When the 
object is a black hole, dissipation is intrinsically built 
in the problem under the form of a one-way membrane: 
the event horizon. Black hole superradiance was in fact 
shown to exist and to trigger interesting phenomena. 
These include black hole bombs, once the black hole is 
surrounded by some confiniM medium, giving rise to ex¬ 
ponentially growing modes [J-S] • The scattering of mas¬ 
sive fields belongs to this category: the mass term effec¬ 
tively confines the field giving rise to floating orbits @ 
or to superradiant instabilities, which extract rotational 
energy away from the black hole (lol - [T^ . The instabil¬ 
ity leads to a depletion of mass and angular momen¬ 
tum from the black hole, resulting in observationally- 
interesting “holes” in the Regge-plane [iM3- As a con¬ 
sequence, the observation of supermassive black holes 
alone can impose stringent bounds on the mass of ul¬ 
tralight bosons dill. Finally, in multi-scalar theo¬ 
ries or theories with a complex scalar, superradiance is 
the responsible for new hairy black hole solutions [13- 
[Hjl (whether or not these solutions are generically ap¬ 
proached as the end-state of the superradiant-instability 
is discussed in Ref. 0). For a Review, we refer the 
reader to Ref. [3|. 

Given the rich phenomenology of black hole superradi¬ 
ance, one might ask whether those results carry through 
to other self-gravitating objects. Ideal stars do not dis¬ 
play superradiance and none of the associated instabili¬ 
ties (except for possible ergoregions instabilities [3]). But 
realistic stars do dissipate energy and they have several 
channels available for it. The original arguments by Zel- 
dovich already made clear that any classical system with 
dissipation will be prone to superradiant scattering QSl 
(see also Ref. [1^ who studied the correspondence be¬ 
tween superradiance and tidal friction on viscous Newto¬ 


nian anisotropic stars). In fact, stars in full General Rel¬ 
ativity were recently shown to disp lay superradiance, as 
long as dissipation is included [^. The purpose of this 
work is to compute explicitly the amplihcation factors 
when a toy model for dissipation is used, and to show 
that, within this toy model, stars are superradiantly- 
unstable against massive scalar-field fluctuations. 

Our results make possible the existence of “hairy” 
stars0, supported solely by superradiance. For dark mat¬ 
ter models comprised of light degrees of freedom, our re¬ 
sults open the door to constrain the interaction cross sec¬ 
tion with ordinary matter or the self-annihilation cross- 
section. 


II. SETUP 

We are interested in describing ultra-light bosonic de¬ 
grees of freedom, described by a scalar minimally cou¬ 
pled to gravity, and interacting with a star. Dissipation 
can be modelled in several - and complex - ways. For 
instance, one can let the scalar couple to spacetime de¬ 
grees of freedom, or one can consider the coupling to mi¬ 
croscopic degrees of freedom, describing the interaction 
of self-annihilation of - for example - dark matter with 
a star. The description of a dissipative system within 
full General Relativity for systems (such as stars) with 
man y in ternal degrees of freedom is a very complex prob¬ 
lem j26l | . 

Instead, we follow Zeldovich’s seminal work and con¬ 
sider a toy model for absorption, by modifying the Klein- 
Gordon equation inside the star and in a frame co¬ 
rotating with the star as 

□ = ( 1 ) 


^ by which we mean new stationary configurations describing a 
star surrounded by a scalar-field condensate 
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Here, is the mass of the scalar field. The a term is 
added to break Lorentz invariance, and describes absorp¬ 
tion on a timescale r ~ 1/a, in a frame co-rotating with 
the star. The constant 1/a has dimensions of time, or 
length in geometrical units. Contact with particle physics 
can be made by considering the particle associated with 
the scalar field and identifying 1/a with the mean free 
path i of the particle inside the star. In other words, for a 
particle with cross section cr to decay into other particles 
via interaction with say a neutron of mass mjv, then 


a = ncr, (2) 

where n is the nucleon number density in the star. One 
can also express the dissipation parameter a with the 
help of the star’s total mass M, in the dimensionless com¬ 
bination 
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Thus, the dimensionless combination can be of order one 
for interesting dark matter candidates [27l - [^ (and ref¬ 
erences therein). The constant a could also be used sim¬ 
ilarly to describe dark matter self-annihilation. 

It is easy to see, following Zeldovich, that if the fre¬ 
quency in the accelerated frame is uj and the field behaves 
as then in the inertial frame the azimuthal 

coordinate is Lp = p' — Qt, and hence the frequency is 
Lo' = UJ — toO. In other words, the effective damping pa¬ 
rameter auj' becomes negative in the superradiant regime 
and the medium amplifies - rather than absorbing- radi¬ 
ation [I|, [ 2 ^ . This procedure can be used to derive the 

wave equation of a scalar field in a rotating black hole 
background (up to first order in rotation), starting from 
that of a scalar field in a Schwarzschild geometry. As we 
will shortly see, it also allows to predict the amplification 
factors of scalar fields in a Kerr background. Thus, while 
only a phenomenological approach to the problem, the 
inclusion of the Lorentz-violating term is well motivated 
and consistent with known results. In Appendix we 
work out another model for dissipation, inspired in the 
theory behind dissipation in fluids [1^. The results are 
perfectly consistent with the toy model (HD adopted in 
the main body of this work. 

As is clear from the previous discussion, superradi¬ 
ance requires only rotation and dissipation. The par¬ 
ticular background geometry or material is only relevant 
to quantitatively describe the phenomena. We will focus 
exclusively on backgrounds describing a constant density, 
perfect fluid star in General Relativity as a proxy for the 
generic case: 


with 

3 \3y/l-2M/R- y/l-2Mr^/R^ J 

e'^’ = ^ v^l - 2M/R - ^\/l - 2MrVi?3 . (5) 

Here, M, R are the star’s mass and radius, and are related 
to the density p via M = ^pR^. 

Notice that we take an exact general-relativistic solu¬ 
tion describing a non-rotating star. Although not self- 
consistent, the effects of General Relativity are subdomi¬ 
nant. The leading effect of rotation is taken into account 
when transforming to a frame at rest with the star, and 
using (P). 


III. SUPERRADIANT AMPLIFICATION 
FACTORS OF MASSLESS FIELDS 



FIG. 1. Superradiant amplification of a scalar field by a rotat¬ 
ing star, where dissipation is modeled through equation ©. 
Here, Ma = 0.1, Mfl = 0.01. We find that the amplification 
factor scales with the dissipation parameter a for small Ma. 
For non-compact, slowly-rotating stars, the amplification fac¬ 
tor (solid lines) is well described by the analytic result ([S|). 
shown as a dashed curve in the plot. 

To solve the Klein-Gordon equation (ED, we use the 
following ansatz for the field, 

$(t,r-,d,</p) = ^^e-“‘+*’"‘^P/(cosi?), (6) 

r 

where Pi are Legendre functions. To solve the equation 
inside the star, we first transform to a frame co-rotating 
with the star (which amounts to replacing ip — >■ ip — fit, or 
UJ ^ UJ — mfl), impose regularity at the center and then 
match the radial function and its derivative. For now we 
focus on zero field mass, ^ = 0. At large distances, the 
solution takes the form 

d/ = Aoute“" + , (7) 


= —e^'^dP -\- 


dr^ 


1 — 2m{r)jr 


r^dfl^ 


where the subscripts denote an out-going and an in-going 
piece (due to the time dependence of the ansatz for the 
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field). Thus, the asymptotic solution allows one to frame 
this in a context of a scattering experiment. We define 
the amplification factor as |^outP/|AnP — 1- Positive 
amplification factors correspond to superradiant ampli¬ 
fication from the star. Negative factors correspond to 
overall absorption by the star. 

Our numerical results are summarized in Fig. [T] for 
Ma = 0.1, Mfl = 0.01. Superradiant amplification does 
exist, as was shown generically in Ref. |25| . and the am¬ 
plification can be significant. In fact, we find that it 
scales linearly with Ma for small Ma and it increases 
significantly as the star’s surface velocity increases. Am¬ 
plification factors can be of order of those around rotating 
Kerr black holes or higher. 

For non-relativistic, Newtonian configurations 
{M/R, Hi? <C 1) the wave equation can be solved 
analytically inside and outside the star in terms of 
Bessel functions. We find a simple analytic expression 
for the amplification factor, 

lAoutP 4aR{Rn-u;R){ojRf+^ 

|Ai„|2 (2; + l)!!(2Z + 3)!! ‘ 

As can be seen from Fig. [T] this relation gives a very 
good approximation to the numerical results for a good 
fraction of the parameter space, including relatively com¬ 
pact stars. This relation is also interesting, as it allows 
one to predict the amplification factor for rotating black 
holes. For black holes, R ~ 2M and 1/a = M is the 
only possible timescale in the problem. For example, the 
above relation predicts that slowly rotating black holes 
in General Relativity amplify I = 1 scalar fields with 

\^^-l = ^Min-ca){2Mcaf . (9) 

On the other hand, a matched-asymptotic expansion cal¬ 
culation in full General Relativity yields approximately 
the same result (the coefficient turns out to be 2/9 ~ 0.22 
instead of 16/45 - 0.355 @,[311). 

IV. SUPERRADIANT INSTABILITY AGAINST 
MASSIVE FIELDS 

As soon as the field is allowed to be massive, superradi¬ 
ance can trigger new phenomena, in particular instabili¬ 
ties: the mass confines the field within a distance of 
the star, whereas superradiance amplifies it. This effect 
was shown to occur for black holes [1 [iM3 , we now 
show that it occurs generically. 

For massive fields, the asymptotic behavior at spatial 
infinity is of the form ^ Thus, it is possible 

that some eigenvalues cj = ujr -|- iuji exist for which the 
field is regular at infinity (for ujji < fj.) and for which the 
mode is unstable (for w/ > 0, cf. Eq. (|6l)). We have used 
a direct integration approach to search for these modes. 

Figure [2] summarizes our results for the instability of 
massive scalars coupled to dissipative stars. For a fixed 


star radius and mass (i?, M) and dissipation parameter 
a, there are unstable modes. At very small values of M^, 
the instability timescale decreases when the scalar mass 
/r increases. We hnd that at low masses /r our results 
are well described by wj oc (w/j — mfl){M with rj ^ 
9 ± 2. However, the real part of the mode ujr ^ /r also 
increases. Thus, for /j, ^ H the m = 1 mode ceases to be 
unstable. This is clearly seen in Fig. [21 which shows that 
the transition from instability to stability coincides with 
a change in sign of w/j — mQ. 


V. DISCUSSION 

Neutron stars or other highly compact stars with dis¬ 
sipative channels may be interesting probes of physics 
beyond the Standard Model, specially of bosonic fields 
with a Compton wavelength of the order of the star size. 
Our working model is extremely crude and more sophisti¬ 
cated ways of including dissipation should be sought for. 
Nevertheless, our work certainly shows that there may be 
a payoff to these efforts: highly spinning compact stars 
develop instabilities a gain st massive fields which may be 
used (see e.g. Refs. @,[l3) to further constrain ultralight 
bosonic fields, while the bosonic clouds formed during the 
process might emit very distinctive gravitational-wave 
signals within the optimal sensitivity band of ground- 
based detectors 0- It is even possible, but it remains 
to be proven in a concrete theory with dissipation, that 
new star-like configurations exist, which would be the 
counterpart to hairy black hole solutions recently studied 
in theories with a minimally coupled massive scalar [13- 
0111 . 
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FIG. 2. Superradiant instability of a massive scalar field coupled to a rotating star, where dissipation is modeled through 
equation ©• Here, R/M = 4, MQ, = 0.2, Ma = 20 (left panel) and M/r = 0.17 (right panel). The transition from instability 
oj/ > 0 to stability loj < 0 coincides with a change in sign of ur — mQ. We find that Mlui scales linearly with the dissipation 
parameter a for small Ma. 


Appendix A: Superradiant amplification factors in a 
viscons gravity analogue 



FIG. 3. Superradiant amplification of a scalar field by a rotat¬ 
ing star, where dissipation is modeled through equation (EH. 
Here, = 0.002, RQ. = 0.5. 

A different way to model dissipation in a gravitational 
system is to use the methods of analogue gravity. For 
example, linear perturbations in an inviscid, irrotational 
flow propagate as fields on a curved space-time. Thus, 
one natural option to model dissipative phenomena is 
through the inclusion of a kinematic viscosity. For a fluid 
at rest, the wave equation for the perturbations is given 
by 

dU = y'"(t>+\^dtV^cl>. (Al) 

Thus, we will here briefly describe how our results change 


when the dissipation inside the star is modeled through 
(EH). We are here concerned only with rotational su¬ 
perradiance in flat spacetime, and we therefore take 
to be the flat-space Laplacian in spherical coordinates. 
The process of transforming to a co-rotating frame and 
the imposition of regularity at the center are identical to 
that of the model studied in the main text. The amplifi¬ 
cation factors were numerically obtained and are shown 
in Fig. [3] Note that for low frequencies it is possible 
to map this analogue model into the previous toy model. 
Equation EH) can be written as 

^1 — = 0, (A2) 

which in the limit w <C 1 is the same as Eq. ([T]), for 
a = This map between a and v allows us to 

obtain the non-relativistic analytical expression for the 
amplification factors from equation ([8|: 

|Ai„|2 3(2Z + l)!!(2/ + 3)!! ' ^ ’ 

It is worth noting that, if one were to try the naive 
black hole limit R = 2M with i/ of order M, the above 
result would never recover the correct behavior for the 
amplification of waves by rotating black holes. The rea¬ 
son, it seems, is that dissipation in non-rotating black 
holes can be best described through o rather than this 
viscous analogue. This is possibly due to the fact that 
the dispersion relation of Eq. ED in the eikonal limit 
leads to an energy-dependent dissipative term [s^. 
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